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Ice accretion codes depend on models of roughness parameters to account for the 
enhanced heat transfer during the ice accretion process. While mitigating supercooled large 
droplet (SLD or Appendix O) icing is a significant concern for manufacturers seeking future 
vehicle certification due to the pending regulation, historical ice roughness studies have been 
performed using Appendix C icing clouds which exhibit mean volumetric diameters (MVD) 
much smaller than SLD clouds. Further, the historical studies of roughness focused on 
extracting parametric representations of ice roughness using multiple images of roughness 
elements. In this study, the ice roughness developed on a 21-in. NACA 0012 at 0° angle of 
attack exposed to short duration SLD icing events was measured in the Icing Research 
Tunnel at the NASA Glenn Research Center. The MVD’s used in the study ranged from 100 
jam to 200 jam, in a 67 m/s flow, with liquid water contents of either 0.6 gm/m 3 or 0.75 gm/m 3 . 
The ice surfaces were measured using a Romer Absolute Arm laser scanning system. The 
roughness associated with each surface point cloud was measured using the two-dimensional 
self-organizing map approach developed by McClain and Kreeger (2013) resulting in 
statistical descriptions of the ice roughness. 


Nomenclature 


Ac 

= 

accumulation parameter 

AOA 

= 

angle of attack 

b 

= 

codebook vectors 

h(i,j ) 

= 

neighborhood function of i to j codebook vectors 

j 

= 

codebook vector index 

LWC 

= 

liquid water content [gm/m 3 ] 

M 

= 

number of codebook vectors 

MVD 

= 

median volumetric diameter [pm] 

N 

= 

airfoil or mean ice shape surface normal coordinate direction 

R 2 

= 

coefficient of determination (regression) 

R d 

= 

high-dimensional data space 

RMH 

= 

roughness maximum height 

** 

= 

the root-mean-square or “standard deviation” roughness height 

r a 

= 

leading edge radius of curvature 

SOM 

= 

Self-Organizing Map 
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= airfoil or mean ice shape surface tangential coordinate direction 
= Standard error of the estimate for regression 
= element of data set 

= local direction angle of manifold through a codebook vector 
= manifold 

= direction angle of surface point relative to manifold direction through winning codebook vector 
= scaling parameter governing neighborhood size 
- learning rate 

I. Introduction 

T HE effects of roughness on ice shapes are well documented. Ice roughness increases the skin friction drag, 
increases the rate of heat transfer from the freezing water layer to the airstream, and induces higher levels of 
turbulence in the airflow than would be observed on mean representative ice shapes without roughness. Further, ice 
roughness increases the droplet capture efficiency of the local surface relative to a smooth surface. 

Capturing the enhanced convective heat transfer caused by the ice roughness is a critical aspect in modeling the 
accretion process. Since the freezing process is exothermic, energy must be transferred from the freezing front to 
the airstream. The local rate at which energy is transferred governs the rate of ice accumulation. Consequently, the 
regions with the highest convection coefficients are expected to exhibit the highest rates of ice shape growth. 

While ice roughness and its geometric properties have been recognized as critical to the ice accretion process and 
to the reduction in iced aircraft performance in icing conditions, describing the nature of ice accretion roughness is 
not trivial. Classical studies of ice roughness during short-duration icing events include Anderson and Shin [1], 
Anderson et al. [2], and Shin [3]. For the afore mentioned measurements, images were acquired of the ice surface 
from multiple directions. Image analysis techniques were then used to determine the dimensions of each roughness 
element and the spacing to its nearest neighboring roughness element. The methods used in the classical studies 
were extremely laborious in terms of manual identification of each roughness element, use of the image analysis 
equipment, and identification of the roughness element in multiple images. Further, the results of the 
characterizations were morphological or deterministic, relating to individual bead height and width, as opposed to 
the statistical representations commonly used to describe randomly rough surfaces. 

The historical investigations of ice roughness also focused on what is known as Appendix C roughness. 
Appendix C icing represents stratiform and cumuliform icing conditions that exhibit supercooled or partially frozen 
droplets with an MVD of 1-50 pm [4]. While Appendix C conditions were the primary focus of ice protection 
system certification in the past, supercooled large droplet or “Appendix O” conditions have become a significant 
issue regarding future vehicle and ice protection system certification [5]. While Appendix O conditions have MVDs 
much larger than Appendix C conditions, 50-1000 pm, the trends for ice roughness development from the historical 
studies of Appendix C conditions are often extrapolated to Appendix O conditions when using ice accretion codes 
such as FEWICE [6]. 

In addition to changing certification requirements, recent improvements in laser scanning approaches have the 
potential to revolutionize the process of ice accretion shape characterization and surface roughness evaluation. Fee 
et al. [7] recently demonstrated measurements of different ice shapes formed in the NASA Icing Research Tunnel 
made using a Romer Absolute Arm laser scanning system. While laser scanning represents a tremendous leap in ice 
shape and roughness measurement ability, laser scans of ice shapes or castings of ice shapes present challenges. 
First, the scans typically result in point densities of 5000 points per cm 2 of scanned surface area. For the leading 
edges of two-dimensional wings (airfoils) used in the NASA icing tunnel, which have 1 8 to 2 1 in. chord lengths, ice 
shape scans with a 10-cm spanwise width typically result in more than 1 million measurement points. Secondly, 
on two-dimensional airfoils, horn structures and the geometrically closed-loop nature of airfoils lead to multi-valued 
functions describing the airfoil surface. The multi-valued function description of airfoils with complex horn shapes 
results in situations where standard regression cannot be employed to describe the entire airfoil at once. To 
evaluate the roughness along an iced airfoil, the mean ice-shape or ice form must be extracted. Thirdly, the radii of 
curvature of horn structures may be on the same order of magnitude as the roughness heights for long-accretion time 
ice shapes. The comparable sizes of horn structures and roughness creates difficulties in discerning the difference 
between ice form changes and roughness features. Finally, most ice shapes, even on straight wings, have 
considerable three-dimensional variations along the span of the wing. 

Clustering methods for non-linear manifold detection enable a rigorous evaluation of ice shape roughness. 
McClain et al. [8] applied the self-organizing map (SOM) approach to characterize the nature of airfoils with rime 
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ice, glaze ice, and multi-horn ice shapes. The primary objective of McClain et al. [8] was to demonstrate how the 
SOM could be used not only to capture the mean ice form or shape, but to expand the SOM approach using multi- 
dimensional statistics to characterize the ice-surface coverage limits. McClain and Kreeger [9] revisited the SOM 
approach and demonstrated a method to unwrap the point cloud using SOM representations and use the distribution 
of the surface points closest to each SOM codebook vector to evaluate the local surface roughness. 

The objective of this effort was to characterize SLD ice roughness in short duration icing events using the Romer 
laser scanner and employing the approach of McClain and Kreeger [9] to create statistical descriptions of the ice 
roughness. In all of the cases, the icing events were ended before enough ice was accumulated on the airfoil 
surfaces to substantially alter the aerodynamic performance. The primary aims of the effort were 1) to characterize 
how SLD roughness develops with increasing accumulation parameter (indicating longer accretion times) given 
fixed cloud properties and 2) to characterize the influence of the MVD on ice roughness provided equivalent 
accumulation parameters. 


II. Experimental Method 

All of the experimental measurements were performed in the Icing Research Tunnel (IRT) at NASA Glenn 
Research Center in Cleveland, OH. A straight, 21-in. (53.3-cm) NACA 0012 airfoil, which spanned the entire 72-in. 
(182.9 cm) test section, was used for the study. The test airfoil was set to 0° AOA for all cases. 

With the exception of the MVD values, the conditions used for the study were based on the conditions employed 
in the historical literature for the Appendix C cases [l]-[3]. All cases were run with a nominal freestream velocity 
of 129.7 knots (66.723 m/s). Table 1 presents the freestream total temperature (T tota i), the frestream velocity, the 
median volumetric diameter, the liquid water content, the accumulation parameter and the freezing fraction. The 
accumulation parameter is a derived quantity calculated using 

LWC -V- At, 

A c = - 

^ ' r a Pice 

Where r a is the leading-edge radius of curvature of the NACA 0012 airfoil, and p ice is the density of ice at the 
freestream temperature. 


Table 1. Cases Used for the SLD Roughness Characterizations 


Run Number 

T total 

(°F) 

V 

(knots) 

MVD 

(H m ) 

LWC 

(gm/m 3 ) 

Af s 

(s) 

Accumulation 

Parameter 

n 

Freezing 

Fraction 

112912.04 

27.8 

129.7 

100.0 

0.60 

27 

0.0703 

0.217 

112912.03 

27.8 

129.7 

100.0 

0.60 

43 

0.1120 

0.217 

112912.02 

27.8 

129.7 

100.0 

0.60 

55 

0.1432 

0.217 

112912.01 

27.8 

129.7 

100.0 

0.60 

75 

0.1953 

0.217 

112812.07 

27.8 

129.7 

100.0 

0.60 

94 

0.2448 

0.217 

113012.06 

27.7 

129.7 

150.0 

0.60 

27 

0.0703 

0.217 

113012.05 

27.7 

129.7 

150.0 

0.60 

43 

0.1120 

0.217 

113012.04 

27.7 

129.7 

150.0 

0.60 

75 

0.1953 

0.217 

112912.06 

27.7 

129.7 

150.0 

0.60 

94 

0.2448 

0.217 

112912.08 

26.5 

129.7 

200.0 

0.75 

19 

0.0618 

0.217 

113012.01 

26.5 

129.7 

200.0 

0.75 

27 

0.0878 

0.217 

112912.09 

26.5 

129.7 

200.0 

0.75 

43 

0.1398 

0.217 

112912.10 

26.5 

129.7 

200.0 

0.75 

75 

0.2439 

0.217 

112912.07 

26.5 

129.7 

200.0 

0.75 

94 

0.3057 

0.217 


Table 1 shows that for the 200-pm cases, a different liquid water content was used than for the 100-pm and 150- 
pm cases. It was desired to use an LWC of 0.6 gm/m 3 for all cases. However, a 200-pm cloud with an LWC of 0.6 
gm/m 3 case was outside the calibration envelope of the spray-bar system. Consequently, the 200- pm cases were run 
at 0.75 gm/m 3 and the spray times were shortened to match accumulation parameters. 
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For each test, the airspeed and freestream total temperature were set and the spray bar pressures were selected to 
provide the appropriate LWC, MVD, and freezing fraction. A thermocouple embedded in the NACA 0012 airfoil 
was used to determine when the airfoil had reached thermal equilibrium with the flow. Once the thermocouple 
reported the static temperature of the flow indicating that it was in thermal equilibrium with the flow, the spray bars 
were actuated and closed after the predetermined spray time. 

Following the completion of the icing spray, the wind tunnel velocity was reduced to approximately 10 knots 
while keeping the static temperature around -4 °C to avoid thawing of the ice shape. The iced airfoil was painted 
using an alcohol-based titanium dioxide paint. The Romer Absolute Arm was placed in the IRT test section 
upstream of the airfoil. A scan was then made of the leading 120-150 mm of the airfoil leading edge. The scans 
were approximately 100-cm wide (in the spanwise direction) and were performed at the center-span location of the 
airfoil. Figure 1 presents a point-cloud acquired during the tests represented as a water-tight surface stereo- 
lithography (.stl) file. 



Figure 1. Point Cloud and Water-Tight Surface Representation of the 113012.4 Case 

III. Ice Roughness Characterization 

The ice shape point clouds were processed using the Self-Organizing Map (SOM) approach developed by 
McClain and Kreeger [9]. The self-organizing map, or sometimes referred to as a Kohonen Map, is a clustering 
method for the detection of non-linear manifolds, which may be curves or a surfaces, in multi-dimensional space 
[10]. SOMs depend on the use of codebook vectors, b, which may also be called codebook points or neurons, to 
represent clumps of data. A set of codebook vectors representing clumps of data points is depicted in on the left side 
of Figure 2. In its simplest essence, self-organizing maps are employed to capture trends of large data sets by 
representing those large data set by a relatively small set of codebook vectors. 

To capture the trends of a large data set, an SOM algorithm begins by distributing a number of codebook vectors 
randomly through the space contained by the data set. The codebook vectors are then sequentially moved in the 
direction of the clump of data points that are closest to the codebook vector. Like most neural network approaches, 
the SOM requires a learning or training process. Over iterative moves, the codebook vectors spread out and settle 
into their local clumps. When a series of codebook vectors are connected in a sequence, as shown on the right side 
of Figure 2, the representative manifold, or curve as is the case in Figure 2, guiding the data through the data space 
is represented by the path or daisy-chain through the codebook vectors. For more detailed information on self- 
organizing maps and their application in ice shape description and roughness evaluation, please consult Refs. [8] and 

[9]. 

When applied to an iced airfoil or wing without sweep or significant spanwise shape changes, the SOM is 
expected to identify a curve in the Chord-Chord Normal plane (or x-y plane when the spanwise axis is placed in the 
z-direction), which represents the mean shape of the iced airfoil. To demonstrate the SOM characterization of an 
iced airfoil, Figure 3 presents an example point cloud and SOM representation used in the study. The surface 
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presented in Figure 3 is the 112912.07 case which exhibited the largest roughness and the highest accumulation 
parameter (0.3057) measured in this study. Figure 3 demonstrates that the blue codebook vectors capture the shape 
or curve of iced surface in the projected x-y space. 



Figure 2 . Self-Organizing Map Representation of Point Cloud [9] 
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Figure 3. Point Cloud and SOM Representation of 112912.07 Case (MVD = 200 um , LWC = 0. 75 gm/m 3 , and At s 

= 94 sec.) which Exhibited the Largest Roughness Values 

Figure 3 also demonstrates that different densities of codebook vectors were used in different regions of the 
airfoil to capture the surface. Following the method presented by McClain and Kreeger [9], a separate SOM 
characterization with a spacing between 1 mm and 2 mm was used to represent the smooth section near the leading 
edge of the iced airfoil. Further away from the leading edge, the spacing between the codebook vectors was 
between 6 mm and 10 mm in arc length along the surface. Further, Figure 3 demonstrates that while a significant 
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amount of ice and roughness was allowed to accumulate on the NACA 0012 airfoil, the ice did not significantly 
change the apparent aerodynamic shape of the airfoil. 

The nature of the SOM method and the positioning of the codebook vectors along a “daisy-chain” enable a 
statistical evaluation of iced airfoil surface roughness. Since the “clumps” of points are distributed about the 
codebook vectors, the deviations of the point measurements in the clumps can be used to evaluate the coverage 
statistics and uncertainty of the codebook vector representation. Figure 4 demonstrates how each surface 
measurement is used to determine a deviation from the spline surface through the control points or codebook 
vectors. 

Figure 4 shows a single surface measurement x J , its closest codebook vector b n . The two neighboring codebook 
vectors along the daisy-chain of codebook vectors representing the manifold are also shown. In the approach used 
for this study, the manifold is assumed to be a first-order manifold in two-dimensional space with the characteristic 
that at each codebook vector, the local slope of the manifold is equal to the central finite-difference evaluated using 
the two closest surrounding codebook vectors. The approach used assumes that all deviations from the manifold are 
normal to the manifold. That is, the deviation of a surface measurement normal to the line through the codebook 
vector with the local slope set by the neighboring codebook vectors is considered the “height” of the surface point 
above or below the local manifold. 
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Figure 4. Metrics of Local Point about a Codebook Vector 


To calculate the local height, of any point, x J , relative to the manifold, the two neighboring codebook vectors are 
first used to calculate the direction of the manifold through b n using 


a b a = arctan 


y h n+1 y h nl 


(i) 


The direction of the x J point from its winning codebook vector relative to the line through the codebook vector with 
the direction a is then found using 


v j = arctan 


y* 


-y h ° 


-v -* b " 


~ a b n 


The normal height the x J point from the line through its winning codebook vector is then determined using 


dN x ' = k --v ) 2 + 6 v ) 2 ] 2 ) 


( 2 ) 


( 3 ) 


The normal height of all the points related to an individual codebook vector may then be used to calculate 
statistics such as those commonly used to calculate traditional roughness parameters. For example, the root-mean- 
square roughness height for a rough surface is traditionally described as 


*' = 






( 4 ) 
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Based on the SOM-manifold description used here, the root-mean- square roughness height is calculated at each 
codebook vector as 


R 


q 




N 


7=1 


1 

2 


( 5 ) 


In traditional roughness studies, the 99.9% roughness maximum height (RMH) based on a Gaussian distribution is 
calculated using 3.09 times the root-mean roughness height. 

2 

( 6 ) 


RMH. a =3.09 R a ==3.09 
b q 


i J 

J X 


7=1 


Where J is the number of surface points for which b n is the winning (closest) codebook vector. The RMH is the 
local 99.9%-maximum roughness height and is evaluated at each codebook vector based on the local or 
neighborhood statistics. The RMH will vary along the surface arc length. 

To evaluate the distance along the manifold representing the mean ice shape, a discrete arc-length approach is 
taken. That is, at one end of the daisy chain, the length of the arc is set to zero or a known value from the stagnation 
point on the airfoil. The position of the next codebook vector is then evaluated as the straight-line distance between 
the two codebook vectors as demonstrated in Eq. (7). 

s b » = V' + [(-v - x b“-> ) 2 + Gv - V- 1 ) 2 ] 2 ( y ) 


Once the surface distance coordinate of each codebook vector is determined, the surface distance coordinate of 
each point cloud measurement may be evaluated based on the location of the surface point’s winning codebook 
vector. Revisiting Figure 3, once the angle of the x J point with respect to the surface manifold through its winning 
codebook vector (y) is known, the surface projection along the manifold is found using 

s xi =5 b „ +ds xi =s b „ + [(, x , -* b n) 2 +(>V -v) 2 ] 2 cos (r x j) ( 8 ) 


The straight-line arc length approach of Eq. (7) is not the only option for evaluating the length of the manifold 
between the codebook vectors. An approach that connects the codebook vectors using the lines through the 
codebook vectors with the directions and slopes evaluated using Eq. (1) may be a more appropriate method than that 
presented in Eq. (7). However, the salient features of any unwrapping process are captured by the method 
employing Eqns. (7) and (8). Those critical features being 1) the distance between the codebook vectors along the 
manifold are calculated and then 2) the projected distance of each surface point from its winning codebook vector is 
calculated relative to the directions tangent to and normal to the manifold passing through the winning codebook 
vector. 


IV. Roughness Measurement Results 

The specific aims of the study were to 1) characterize how SLD roughness develops with increasing 
accumulation parameter given fixed cloud properties and 2) characterize the influence of the MVD on ice roughness 
provided equivalent accumulation parameters. In the following sections, each aim is discussed separately. 

Accumulation Parameter (Icing Event Duration) Investigation 

Each of the ice shape scans described in Table 1 consists of more than one million surface points. 
Consequently, reducing the point clouds and understanding the implications of the reduced roughness data is 
inappropriate without an appreciation for the rough surface topographies captured in the study. Figure 5 presents a 
scatter plot of the unwrapped surface topographies for four of the 200-pm cases. Since a two-dimensional SOM 
approach was used to characterize the point clouds, the spanwise (z) coordinate of each point in the original point 
cloud was not modified for a given scan. Approximately 20-mm strips from the spanwise center of the scan were 
extracted for Figure 5, and points only on one side of the airfoil were plotted (i.e. the “0” distance point on the x-axis 
represents the stagnation point of the iced airfoil). 
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The color map for Figure 5 was set based on the elevation values of the 94- sec. case to provide a consistent 
comparison of the accretion time-development of roughness. The elevation changes shown in Figure 5(a) for the 
27-sec. case are barely noticeable given the scaling mandated by the 94-sec. case. Because of its very small 
roughness values and because of the lack of resolution for the 27-sec. case, the 19-sec accretion case at 200 pm and 
0.75 gm/m 3 (Run 112912.08) was omitted from Figure 5. As the accretion times grow and the accumulation 
parameter increases, the surface features become more and more pronounced in Figure 5. 
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Figure 5. Surface Topography Visualization along the Airfoil Surface for the 200-um, 0. 75 gm/m 3 Cases 



Figure 6, 7, and 8 present the maximum roughness height (RMH) measurements, calculated as three times the 
root-mean- square roughness height, along the surfaces for the 100-pm, 150-pm, and 200-pm cases, respectively. In 
Figures 6, 7, and 8, the “zero” mean shape surface distance represents the location of the stagnation point on the iced 
airfoil. Positive values of mean shape surface distance represent surface distance variations along the right side of 
the airfoil when an observer is aligned with the spanwise axis of the airfoil and is looking into the direction of the 
oncoming flow. Negative values represent surface distance along the left side of the airfoil. Figures 6, 7, and 8 all 
show that the roughness values along the leading edge-region become more and more rough with increasing ice 
accretion time or accumulation parameter. 

However, Figures 6, 7, and 8 present some issues in quantifying the variation in ice roughness with accretion 
time. First, there is not a defined region where the RMH values are constant in the direction of the flow. The 
magnitude of the RMH varies continuously along the surface. Secondly, Figures 6-8 demonstrate that the roughness 
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variations are not symmetric even though all of the cases were run for a symmetric airfoil at 0° AO A. The value of 
roughness at a location on the right ice surface generally does not equal the roughness at the reflected location on the 
left ice surface of the airfoil. Finally, in investigating the roughness variations with accretion time, the location of 
the roughness comparison matters. For example in Figure 6, the surface maximum RMH exhibited by the 100-pm, 
75-sec. case is actually larger than the 100-pm, 94-sec. case. However, at ±60 mm from the leading edge, roughness 
values appear to progress proportionally to the ice accretion time. 



Mean Shape Surface Distance (mm) 

Figure 6. Roughness Measurements for MVD = 100 jum, LWC = 0.6 gm/m 3 



Mean Shape Surface Distance (mm) 

Figure 7. Roughness Measurements for MVD = 150 jam, LWC = 0.6 gm/m 3 



Mean Shape Surface Distance (mm) 

Figure 8. Roughness Measurements for MVD = 200 jam, LWC = 0. 75 gm/m 3 


Given the asymmetry and surface roughness variations for each case, the SLD cases were compared in their 
development versus accretion parameter using the average-maximum RMH. That is, for a given experimental test 
case, the maximum local roughness height on the right side was averaged with the maximum local roughness height 
on the left side of the airfoil. 


RMH 


max 


RMH max Right + RMH mSiX Left 
2 


( 9 ) 


Figure 9 presents the airfoil RMH max values as a function of accumulation parameter for the SLD cases. Figure 9 
demonstrates a linear relationship between the airfoil maximum roughness height and the accumulation parameter. 
Interestingly, only one data point is an outlier with a deviation of 0.4 mm from the predicted maximum RMH, while 
two of the 100 um cases (the A c = 0.0703 and A c = 0.1 120 cases) are not visible in Figure 9 because they are covered 
by blue circles of the 150 micron cases. 
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Figure 9. Airfoil Maximum Roughness Height as a Function of Accumulation Parameter 


The results from Figure 9 differ from the findings of Shin [3] who studied Appendix “C” roughness on a 53.3 
cm (21 in.) NACA 0012. Shin [3] exposed the NACA 0012 to clouds with MVD = 20 pm and LWC = 0.5 gm/m 3 
and with velocities of 67.1 m/s or 111.8 m/s. Images taken looking down the spanwise axis at the surface of the 
airfoil were then used to measure the roughness element heights. The data reported in Shin [3] are included in 
Figure 9 and demonstrate that while an increase in roughness levels occurred at low accumulation parameters, the 
“Appendix C” roughness measurements plateaued at 0.6 mm. Further, the roughness heights increased less quickly 
in the study of Shin [3] than in the current SLD study. 

Many factors are likely contributors to the difference in the rates of roughness growth and the limits to 
roughness growth noted by Shin [3] and those found in this study. The most significant difference between the 
measurement of Shin [3] and those of this study is the significant difference in the MVD’s explored. Further 
research using laser scanning of “Appendix C” ice roughness is required to determine whether or not the MVD’s 
employed caused the difference between the ice roughness behaviors. Other suspects for the difference between the 
measured roughness growth rates include 1) the difficulties in performing image analysis when observing small 
surface features from long distances using a zoom lens, 2) the arbitrary nature of picking out “representative” 
roughness elements from a field, and 3) the lack of detail in the Shin [3] study of where, in terms of distance from 
the stagnation point, the roughness element was measured. 


MVD Comparison 

In Figure 9, two accumulation parameters exist with measurements made for each of the three MVD values: A c 
« 0.244 and A c « 0.070. Studying the cases at these two accumulation parameters provide a deeper investigation 
into the influence of the MVD on ice roughness statistics and topography. 

Figure 10 presents the local RMH variations along the airfoil surfaces for each of the three MVD cases for the 
lower accumulation parameter 0.070). Figure 10 shows that there are differences between the three cases at the 
lower accumulation parameter. That is, there is some separation between the three cases especially on the left-side 
of the airfoil. However, when the average values between the left and right side are compared, especially at the 
maximum values presented in Figure 9, the differences between the cases becomes negligible. 

Figure 1 1 presents the local RMH variations along the airfoil surfaces for the higher accumulation parameter (« 
0.244). Figure 1 1 along with Figure 9 show that when comparing the surface maximum RMH values, the maximum 
RMH values are essentially independent of the MVD. However, Figure 1 1 suggests that downstream of the location 
of the surface maximum RMH values, variations in surface roughness statistics are affected by MVD. In Figure 1 1, 
the 100-pm case is generally less rough than the other cases downstream of the point of local maximum RMH 
values. While a quantitative assessment is not provided, the smooth region also appears larger for the 100-pm case 
than for the other two cases. 
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Figure 10. Roughness Comparison for an Accumulation Parameter of Approximately 0.070 (Shorter Event) 



Figure 11. Roughness Comparison for an Accumulation Parameter of Approximately 0.244 (Longer Event) 


Figure 12 presents the topographies of the 0.244 accumulation parameter cases. While the maximum RMH 
values are very similar, Figure 12 reiterates the findings from Figure 11: for the 100-pm case, the roughness starts 
further downstream and decays faster than for the other two cases. 

V. Conclusions 

A series of experiments was performed to quantify the development of roughness on iced airfoils during 
supercooled large droplet (SLD) icing. The measurements were performed on a straight NACA 0012 at 0° AO A 
and at conditions matching historical studies of “Appendix C” icing roughness. The resulting ice shapes were 
scanned and analyzed using the method of McClain and Kreeger [9] for roughness evaluation. The primary findings 
of the study were: 

1. Ice roughness varies continuously along the iced surface of the airfoil. Because of the continuously varying 
properties, comparing ice accretion roughness must be performed using the airfoil maximum values or at a 
specific distance along the iced airfoil surface. 

2. For the three MVD’s chosen for the study, the airfoil maximum ice roughness increased proportionally to the 
accretion time or accumulation parameter. This finding deviates from historical Appendix C studies which 
found that roughness height became capped at around 0.6 mm as the accumulation parameters increased. 

3. While the airfoil maximum roughness was found to vary primarily with accumulation parameter, the cloud 
MVD influenced the roughness variation along the surface of the airfoil. As the MVD decreased and for the 
longer accumulation parameter studied, the roughness region started further from the leading edge and 
decayed faster as the distance from the leading edge increased. 

The results of the study provide significant insights into the development of roughness during ice accretion on 
airframe surfaces in SLD conditions. The new insights also represent a substantial resource for modeling the SLD 
ice accretion process using codes such as LEWICE. 
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Figure 12. Topographical Comparison of Surface Features for Cases with an Accumulation Parameter of 0.244 
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